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ABSTRACT 


Silicon photodiodes have proved to be excellent detectors in the visible wave length range. On the other 
hand, high-performance photodetectors development remains needed. This research paper will focus on 
the calculation and improvement of photogenerated current of PN silicon photodetectors enhanced by 
the impurity photovoltaic effect. Various parameters that affect the behavior of the photogenerated 
current will be analyzed and studied using simulation techniques, namely MATLAB. Other aspects for the 
improvement of the device will include the responsiveness and the extension of the operating 
wavelength range. Different types and densities of impurities will be considered, analyzed and discussed. 
Four metallic impurities are used in the simulation, which are namely boron, indium, gallium, and 
aluminum. On the issue of validation of the data obtained from in this paper a practical investigation was 
conducted in a separate research work which involved laboratory preparation of silicon samples with 
indium or aluminum. Only one figure out of many will be presented in this paper. It will be seen that the 
variation of the photogenerated current with wavelength has a similar profile, which indicates the good 
agreement of the results 

© 2013 Elsevier Ltd. All rights reserved. 
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1. Introduction 


platforms for the telecommunication application in the visible 
range of wavelength remains up to date not complete. Recently, 


Silicon photonics development and progress can be considered 
as the most dynamic in the field of integrated optics, with a vast 
production of commercial products [1]. Silicon photodiodes have 
proved to be very effective and efficient photodetectors in the 
range of visible wavelengths. However, the development and 
production of high-performance photodetectors on silicon CMOS 
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however, research literature shows that some silicon photo detec- 
tors in the near infrared range have been proposed and integrated 
for power monitoring and optical interconnection [2-4]. On the 
other hand, recent rapid progress in the technological processes 
for the use of silicon-on-insulator (SOI) substrates has led the way 
for the generation of a wide range of reliable and efficient fully 
CMOS-compatible optical photonics components. These compo- 
nents include high speed modulators, low-loss waveguide, high-Q 
resonators, couples and optically pumped lasers [5-10]. Silicon 
photodetectors with high speed resonant cavity, using reflecting 
silicon on insulator substrates have been reported in some 
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research works [11]. These detectors seem to have 40% quantum 
efficiency and a surface with 90%refelction at 860nm and 
response time of 29 ps. Bulk photodetectors can be classified as 
the oldest and most thoroughly understood generation of silicon 
optoelectronic components. They are now commercial products 
with operating range of wavelength below 1100 nm with the 
occurrence of band-to-band absorption. However, due to its 
transparency silicon is not considered a suitable material for 
photodiodes incorporated in photonics circuits operating with 
wavelength exceeding 1100 nm. In order to capitalize on the low 
cost of standards Si-CMOS processing technology numerous 
photodetectors have been suggested and developed employing 
different physical effects. These effects include, mid-band gap 
absorption (MBA) [12], surface-state absorption (SSA) [13], internal 
photoemission absorption (IPA) and two-photon absorption (TPA) 
[14]. It is very common that semiconductors suffer from extrinsic 
and intrinsic defects. The extrinsic defect is associated with foreign 
atoms from impurities. Intrinsic defects on the other hand are due 
to crystalline defects. Impurities, in many cases, are intentionally 
introduced as dopant or recombination centers so as to enhance 
the fabrication of semiconductor as devices. Unintentional defects 
are introduced during the phase of crystal growth or device 
fabrication. Defects can be classified into two categories, shallow 
impurities and deep centers. Shallow impurities are characterized 
by levels with binding energies much smaller than that of the 
energy gap of the hosting semiconductors. The category of deep 
centers defects are characterized by electronic level located near 
the middle of the band gap. Defects termed generation-recombi- 
nation centers are deep in the band gap. These defects become 
recombination centers when the semiconductor has carrier den- 
sity lower than its equilibrium value. It must be noted that in the 
case of single crystal semiconductors, deep level impurities are 
commonly generated by metallic impurities. 

The inability of the silicon photodetector absorbing photons 
possessing energies less than the band gap (E,) and the recombi- 
nation life time (s) of the photogenerated carriers, limit the 
detector efficiency. However, deliberate inclusion of impurity 
atoms or defects creates energy states within the band gap which 
helps a significant portion of sub-band gap photons to be 
absorbed. This mechanism is known as the impurity photovoltaic 
(IPV) effect, which can be considered as a promising approach for 
the improvement of photodetector efficiency [15,16]. The energy 
states within the band gap behave as additional recombination 
centers, which, unfortunately, severely lower the carrier life time. 
On the whole, non-mid gap and deep level impurity atoms are 
more suitable as they create some balance between the lifetime 
decay and the access to the sub band gap photons. Utilizing this 
concept, indium (In) can be suggested for the IPV effect as it can 
generate energy levels at around 157 meV above the Si valence 
band [15]. A significant improvement in short circuit current, up to 
5 mA/cm?, can be achieved when using a modified Shockley- 
Reed-Hall recombination model [17]. The assumption with the use 
of the model is infinite mobility and the concentration of indium is 
sufficiently high, around 10'”/cm? so as to compensate the n-type 
base [15]. The validation of the simplified Shockley-Reed-Hall 
statistics used for modeling carrier life time in crystalline silicon 
has been thoroughly investigated [18]. The smallest number of 
recombination center density that ensures that the simplified SRH 
model remains accurate is determined. This enables the lifetime- 
based stroboscopic techniques to safely be used to determine the 
number. Improving the efficiency of photodetectors is investigated 
through combined impact of impurity photovoltaic effect and 
compensation [19]. The IPV effect is produced by a deep level 
impurity, while shallow level impurity is employed to achieve the 
compensation. It must be noted that the deep level impurity must 
be kept moderate by comparison to primary dopant in order to 


prevent any deterioration in the carries lifetime. Overall efficiency 
enhancement of photodetectors can be achieved by the improve- 
ment of the carrier lifetime through compensation. In this research 
paper, the parameters that affect the behavior of the photogener- 
ated current will be analyzed and studied using simulation 
techniques. The simulation will be focus on the calculation and 
improvement of photogenerated current of PN silicon photode- 
tector enhanced by impurity photovoltaic effect. The improvement 
of the device will be mainly, concentrated on the device respon- 
sivity and the extension of the range of the operating wave length. 
Different types and densities of impurity will be considered, 
analyzed and discussed. Four metallic impurities are used in the 
simulation, which are namely boron (B), indium (In), gallium (Ga) 
and aluminum (Al). 


2. Modified Shockley-Read-Hall model 


In a standard Shockley—Read-Hall (SRH) model [20], electrons 
and holes be captured by impurity level and thermally excited 
from the impurity level. The six impurities related transition 
included in the modified SRH model is shown in Fig. 1. Under 
steady state condition the impurity level occupancy does not 
change with time, so that the net rate for electron capture (E;ap) 
by the impurity is balanced by the net rate of holes capture (Heap). 
The net recombination rate via an impurity level can be expressed 
as [21] 


np—ny py 
E Tnolp + pï] + Tpo[N + ni] , 


U(x) (1) 


where n and p are the electrons and holes concentrations and tno 


and tpo are SRH low injection lifetimes for electrons and holes and 
they can be expressed as follows: 
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Fig. 1. Modified Shockley-Read-Hall (SRH) transitions model. 
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Fig. 2. Energy levels of some impurities appear at the energy gap of silicon 
material. 
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(3) 


ny = N1 + TnoSnmax 
*— 
Pi =P1 + Tpo&pmax 


where n; and p4 are the impurity electrons and holes concentra- 
tions and N, is the impurity concentration. C, and C, are the 
electrons and holes capture coefficients and Zjmax is the optical 
emission rate of electrons from the impurity level which is fully 
occupied. Zpmax is the optical emission rate of holes from the 
impurity level which is completely empty. The basic idea of Eq. (1) 
is that for some impurities incorporated into a photodetectors, the 
recombination rate (U) may be positive or negative. 

As far as SRH recombination mechanism is concerned, the 
impurity acts either to increase recombination, i.e., U is positive, 
that lead to decrease the carrier life time, consequently decreasing 
current generation, or to increase current generation in the device 
when U is negative. In order to increase the current generation by 
making U negative, this requires 


np—n} pi <0. (4) 


So, n{p} must be as large as possible. Using Eq. (3), Eq. (4) can be 
rewritten as 


np—(y + tnoSnmax) x (Pı + Tpo& pmax) < 0 (5) 


For midgap impurity levels n,=p,=n;, therefore placing strong 
requirements on both photoemission terms in Eq. (5). 

If impurities are added to the semiconductor material, some 
new energy levels will appear inside the energy gap, as shown in 
Fig. 2 [20]. These new energy levels will be distributed throughout 
the energy gap depending on the type of impurity. Some impu- 
rities produce more than one energy level. Adding such impurities 
will extend the range of absorbed wavelengths by producing sub 
energy gaps inside the original energy gap. Photons of energy less 
than the energy band will be absorbed by these sub energy gaps. 
Electron-hole pair generation occurs in two stages, as shown in 
Fig. 1. One of these two stages determines the rate and type of 
generation. The location of energy level and the stage of genera- 
tion must be taken in account when the impurity is chosen. 

A shallow impurity level gives rise to an asymmetry between 
nı and pj. For the case at hand of shallow acceptor p,>>n,, which 
means that only one of the photo excitation process is crucial for 
nip; to be large, i.e., taoSnmax Should be large. The smallness of 
nı indicates that the electron photoemission term Zpmax iS crucial 
for the recombination rate to be negative, thus 


nÏpř=p: Tno&nmax: (6) 


Znmax and ¥pmax Can be expressed in terms of cross sections of 
optical emission (copt) and photon flux density (ppn) as follows: 


Amax 
Enmax = Ne Í Gnopt Ppn(X; A) da (7) 


‘min 


Amax 
Spmax = Ne | Opopt Pon(X, A) da, (8) 
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Fig. 3. Lambertian model to enhance light trapping. 


where gnopt and gpopt are the cross sections of optical emission for 
electrons and holes respectively. 

The Lambertian model shown in Fig. 3 is used to explain the 
light trapping in the device [18]. 

It is suggested that the surface is formed from inverted 
pyramids to disperse the light all over the detector body. Also, 
the device is coated with antireflection material at the front layer 
and reflection material at the rear layer. The angle of incident is 
60° with the vertical to the device surface. The general expression 
of photon flux density at distant x from the surface is 


1 + e74 W-x)ar Rp 


1-e-4wa R, Rr x e Beal), (9) 


Ppp (X, A= ad —Rext(4)) 
where q; is the total absorbance coefficient, Rext is the external 
reflectance at front layer, R; is the internal reflectance at front layer 
and Rp is the internal reflectance at the rear layer. For typical 
working device of no external reflectance (Rex:=0) and full 
internal reflectance (Rp=1), knowing that Rr=1-1/N?, N is the 
refractive index of silicon which is equal to 3.45. Eq. (9) will be 
reduced to @pn(A)=2N? Gext(A) =25@ext(A). That means the internal 
photon flux is enhanced by the factor 25. 


2.1. Photogenerated current simulation 


The photogenerated current density (Jpn) inside the device due 
to impurity doping is given by [22] 


xj 
Jm=4 | -UOJ de, (10) 
Oo 

where q is the electronic charge and x; is the junction depth. 

Substituting Eq. (1) in Eq. (10) yields 

. np-nīpi 

=- dx 

Jon af Trop + PY) + poln + NY 


Substituting Tno» tpo, N{ and př from Eqs. (2), (3), (7) and (8) in 
Eq. (11) and knowing that 


(11) 


ny =g,Nce-E-E/KT and p) = 7 Nye EED as 
t 


where g, is the impurity level degeneracy, N. and N, are the 
effective density of states in the conduction and valence bands 
respectively, E< and Ey are the conduction and valence band edges, 
E, is the effective impurity level, k is the Boltzmann constant and 
T the temperature in Kelvin (see Eq. (13) below). 

Knowing the values of gi, Ne Nv, N, p, Ch, Cy, on op and 
x;=0.3 um, Eq. (13) can be solved numerically using the MATLAB 
program and the values of the photogenerated current can be 
calculated. The photogenerated current depends on the photon 
flux density, which is a function of spectrum wavelength, impurity 
concentration, and its depth inside the device and the device 
reflectivity. 


3. Results and discussions 


All relevant parameters of the elements, which have been used 
in the study and soluble in silicon material, are summarized in 
Table 1 [21, 23-26]. Different types of impurities with different 
concentrations were used to investigate the effect of these para- 
meters on the performance and behaviors of photodetectors. The 
investigations were carried out on samples doped with impurity 
around 0.3 um deep from the surface using simulation techniques, 
mainly MATLAB Programs. 

The variations of photogenerated current with respect to 
wavelength of boron doped samples without front antireflection 
and without rear reflection are shown in Fig. 4. It can be deduced 
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Table 1 
Impurities parameters soluble in silicon material. 


6- N:= 102 1 
0.2 


Impurity parameter In Ga Al B Reference 
Energy level E; (eV) E,+0.16 E,+0. 67 E,+0.57 E,+0.54 [23] 
Degeneracy gt 4 4 4 4 [24] 
Holes density pı (cm) 1.6 x 101° 6.3 x 10” 8.58 x 10” 1.36 x 108 [21] 
Electrons density nı (cm?) 9.14 x 10? 2.33 x 107 1.71 x 10? 1.1 x 10? [21] 
Hole capture coefficient C, (cm s“') 8x 108 6.67 x 107 2x 107 1.11 x 107 [24] 
Electron capture coefficient C, (cm s~!) 2x10" 1x107" 1x 10716 8.3 x 10-7 [25] 
Hole optical emission cnopt (cm?) 0.5 x 10-6 2.6 x 10718 8x 10716 14x 10718 [24, 25] 
Electron optical emission opop: (cm?) 9x 10-8 10 x 10-9 5.3 x 10-9 3.19 x 10-9 [26] 
Solubility (cm?) 1x 10'8 4x10" 2x 10% 6 x 1070 [26] 
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Fig. 4. Variation of photogenerated current with wavelength of boron doped 
sample without front antireflection and rear reflection. 
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Fig. 5. Variation of photogenerated current with wavelength of indium doped 
sample without front antireflection and rear reflection. 


that increasing the impurity concentration improves the photo- 
generated current and also widen the operating wavelength range 
toward near infrared region. More specifically the current starts to 
rise when the concentration of the impurities exceeds the 
value 10'%/cm? and exhibit peak value at N;=10'°/cm? when the 
concentration of the impurities almost equal to the electron 
concentration of the primary dopant in the window layer. Beyond 
this value of concentration, the photogenerated current tends to 
decline in value. 


xj np—[ny + (1/Cn NNi f7* Zonopt(A)Ppn(X, A) dAHIPi + 1/CpNdANe J" Zo pope(A)Ppn(X, 2) dA] 
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Fig. 6. Variation of photogenerated current with wavelength of gallium doped 
sample without front antireflection and rear reflection. 
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Fig. 7. Variation of photogenerated current with wavelength of aluminum doped 
sample without front antireflection and rear reflection. 


be positive when the values of nf and pł decrease, this conse- 
quently indicates that the process at hand is a recombination 
process rather than a generation process. The sample exhibited 
peak responsivity at 1.05 um wavelength, but it decreases to 
0.4-0.7 um at lower and higher impurity concentrations. 

Figs. 5-7 display the variation of the photogenerated current 
with respect to wavelength without front antireflection and with- 
out rear reflection for different impurity concentrations of In, 


Jpn = -q 


The explanation for this is the decrease of carrier lifetime. On 
the whole, carrier lifetime is inversely proportional with the 
concentration of the carrier, according to Eqs. (2) and (3). U will 


0 1/CrN[p F P14C1/CpNr){Ne fag 2G popt(A)Ppn(X, 4) da} F ad /Cn Nd {n +N + (1/CoNe){Nt p= Zonopt(A) Ppn(X, A) da}} 


(13) 


‘min 


Ga and Al. On the whole, it can be noted that the responsivity will 
exhibit peak values at around 1.05 um for all samples. Indium 
doped samples had the highest photogenerated current at 
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Fig. 8. Variation of photogenerated current with wavelength of boron doped 
sample with full front antireflection and rear reflection. 
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Fig. 9. Variation of photogenerated current with wavelength of indium doped 
sample with full front antireflection and rear reflection. 
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Fig. 10. Variation of photogenerated current with wavelength of gallium doped 
sample with full front antireflection and rear reflection. 
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Fig. 11. Variation of photogenerated current with wavelength of aluminum doped 
sample with full front antireflection and rear reflection. 


10'’/cm? impurity concentration [27,28]. Gallium and alumium 
doped samples had the highest photogenerated current at 
10!8/cm? impurity concentration. All these results are in 
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Fig. 12. Variation of photogenerated current of the laboratory prepared samples: 
(1) indium and (2) aluminum, with wavelength at -3 V reverse voltage. 


agreement with what was reported by Sclar [29]. It can be 
concluded that the value of the photogenerated current is varied 
with the varation of impurity type and its concentration. The 
varation of peak value of photogenerated current with regard to 
the impurity concentration is attributed to the difference of 
impurity solubility as stated in Table 1. Boron has the highest 
solubilityof 6 x 107°/cm?, while indium has the lowest at 101°/cm?. 
Impurity with the lowest solubility tends to have highest cross 
section optical emission (Snop) which pushes the rate of generation 
(Znmax) Up, as can be deduced from Eqs. (7) and (8). This will lead 
to the increase of the photogenerated current [26]. It is abvious 
now that the solubility of the impurities limits the employment of 
high impurity concentration in photodetectors. Indium concentra- 
tion beyond 10'8/cm? causes the photogenerated current to 
decline, while the current declines for gallium and aluminum at 
concentration 10'9/cm? and at 107°/cm? for boron impurities. 

The varaition of the simulated photogenerated current versus 
the wavelength with antireflection layer at the front and reflection 
layer at the rear are taken into account are shown in Figs. 8-11. It is 
clearly evident that the photogenerated current has increased by 
almost ten-fold for all samples. The reason for this is related to the 
enhancementof the internal photo flux density at the expense of 
the external density. It is very interesting to note the extention of 
the oprating wavelength range to 1.1 um at which the photogen- 
erated current reaches the peak. This may support the idea of 
using some kind of light concentrator to push the detection range 
of silicon photodetectors far into infrared. 

The photoelectronic properties of practically fabricated silicon 
photodetector doped with indium or aluminum were investigated 
and compared with simulated results. It can be deduced that high 
photogenerated current were produced by silicon samples doped 
with indium of thickness 1500 A or aluminum of thickness 2000 A. 
These thicknesses correspond to around 10" /cm? for indium and 
around 10!8/cm? for aluminum. 

Fig. 12 shows the variation of photogenerated current of the 
laboratory prepared samples with wavelength at -3 V reverse 
voltage. In general the profiles of the curves are the same as that 
produced by the simulation. It can be noted that the photogener- 
ated current for aluminum samples is less than that for indium 
samples despite the fact that the impurity concentration of 
aluminum is higher. This can be explained as result of the fact 
that optical emission cross section of indium is higher than that for 
aluminum. 


4. Conclusion 


Research in the development of high-performance photodetec- 
tors with sufficient photogenerated current and wide visible 
wavelength suitable for such applications as optical interconnec- 
tion and power monitoring is of great interest nowadays. 
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Simulation techniques are powerful and useful approaches to 
study the enhancement of the performance of PN silicon photo- 
detectors. The results obtained indicate that increasing the impur- 
ity concentration improves the photogenerated current and 
moreover widen the operating wavelength range toward near 
infrared region. The increase in the current exhibits peak value 
that depends on the level of impurity concentration and the type 
of impurity used. However, it can be deduced that, on the whole 
and for all four types of impurities, the responsivity of the device 
reaches peak values at 1.05 um wavelength. Indium doped samples 
produced the highest photogenerated current with the lowest 
impurity concentration of 10'’/cm? as compared with boron, 
gallium and aluminum. These simulated results are for samples 
without front antireflection and without rear reflection. On the 
other hand simulated photogenerated currents versus the wave- 
length for samples with antireflection layer at the front and 
reflection layer at the rear shows that the photogenerated current 
can be increased by almost ten folds for all samples. This is due to 
the enhancement of the internal photo flux density at the expanse 
of the external density. 
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